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SUMMMIY 


A  one- dimensional  model  for  the  Boltzmann  equation  of  kinetic 
theory  of  gases  is  proposed,  in  order  to  study  the  validity  of  the  usual 
approximation  method  and  the  behavior  of  the  solutions  near  the  free  mole¬ 
cule  limit.  The  interaction  term  is  very  similar  to  the  Fokker  Planck 
term.  The  conservation  of  mass,  momentum  and  energy  results  from  the 
particular  form  assumed  for  the  interaction  term  and  the  only  equilibrium 
solutions  are  shown  to  be  of  Maxwell-type.  Hydrodynamic  equations  similar 
to  the  Navier-Stokes-  Fourier  equations  are  derived  in  the  tiisua|:  wa^  under 
the  assumption  of  small  gradients.  Shock  waves  are  shown  to  exist  in 
supersonic  flows. 

The  distribution  function  is  obtained  in  the  case  of  a 
spatially  uniform  condition,  and  Che  existence  of  an  infinite  number  of 
relaxation  modes  and  corresponding  relaxation  times  is  shown. 

The  linearized  versions  of  the  model  equation  are  established 
with  particular  emphasis  on  the  simplest  problems  of  heat  conduction  and 
wave  propagation  and  the  usefulness  of  Fourier  transformation  with  respect 
to  the  velocity  component  is  indicated. 

Finally,  the  behavior  of  the  solutions  near,  the  free  mol¬ 
ecule  limit  is  briefly  investigated  for  steady  state  conditions.  Exist¬ 
ence  of  a  singular  behavior  for  slow  molecules,  analogous  to  boundary 
layer  phenomenon,  is  indicated  and  the  equation  is  shown  to  reduce  to  a 
singular  parabolic  equation  studied  by  Gevrey,  in  the  region  of  interest. 

The  leading  terms  in  the  deviation  from  free  molecule  data 
ate  found  to  be  of  the  order  of  the  cube  root  of  an  interaction  parameter, 
playing  a  role  analogous  to  the  Knudsen  number. 
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1 .  Usefulness  and  Limitations  of  a  Mathematical  Model 

The  mathematical  difficulties  arising  in  the  treatment  of 
the  Boltzmann  equation  for  the  velocity  distribution  function  f  in  a  mon¬ 
atomic  gas  at  moderate  or  low  densities  are  known  to  be  formidable,  a  fact 
which  has  precluded,  so  far,  the  obtaining  of  any  exact  particular  solution 
of  interest  for  conditions  far  removed  from  equilibrium. 


(ref.  1) 


The  Boltzmann  equation  may  be  written  in  the  standard  form 


In  this  equation,  the  distribution  function  f  is  a  (unknown) 
function  of  the  variables  t  (time),  Xj^,  x^,  x^  (space  coordinates) 
(molecular  velocity  components) .  In  the  integral  of  the  second  member, 
the  so-called  collision  integral,  there  appear  the  variables  b  and  which 
define  the  geometry  of  a  particular  collision  (b  being  the  distance  of 
closest  approach  in  the  absence  of  interaction  and  C  an  angle  defining 
the  position,  in  space,  of  the  plane  of  relative  motion  of  the  molecules 
with  respect  to  the  center  of  mass  of  the  system) .  The  dashes  indicate 
the  velocity  components  of  the  partner  molecule  in  the  encounter.  The 

stars  indicate  the  values  of  velocity  components  after  the  encounter,  so 

,  rx  -CA  rr'X  ~Ci^ 

that  >  S  2. ;  S?  ,  S  1  )  ^2.  functions  of  the  initial  velo¬ 

city  components  and  of  b  and  .  Stars  and  dashes  attached  to  the  symbol 
f  indicate  which  velocity  components  should  be  taken  as  arguments.  The 
Integral  is  extended  to  all  possible  encounters  and  all  possible  partners. 


The  velocity  distribution  function  f  is,  of  course,  defined 
in  the  usual  manner  as  indicating  the  nuinber  density  of  molecules  in  con¬ 
figuration  space  (i.e.  the  space  with  coordinates  x^^,  x^,  >^2.  ^  ^3  ) » 
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so  that  the  number  dn  of  molecules  in  an  infinitesimal  volume  dXj^,  dx^, 
dx^  of  space,  centered  at  the  point  (x^^,  x^,  x^)  which  have  molecular 
velocity  components  in  the  ranges  (  ^ 

(  +  is  given  by 

It  is  seen  that  all  the  possible  mathematical  complications 
are  represented  in  eq.  1;  this  equation  being  simultaneously  non  linear, 
partial  differential,  integral  and  functional.  It  is  ntJt  surprising, 
therefore,  that  most  of  the  available  information  concerning  its  solutions 
refer  almost  exclusively  to  two  extreme  limiting  cases:  the  near  equilibrium 
case  and  the  free  molecule  regime  (which  is  somewhat  trivial  in  that  the 
second  member  of  1,  which  contains  all  the  mathematical  difficulties.  Is 
neglected) . 


In  the  near  equilibrium  case,  f  differs  but  little  from  the 
well  known  Maxwellian  velocity  distribution,  given  by  (ref.  1) 


Here,  n  represents  the  number  density  of  molecules  (in  the  ordinary  sensf^ 
i.e.  In  physical  space),  m  the  mass  of  a  molecule,  T  the  absolute  tempe- 
ature,  xl  the  macroscopic,  or  bulk,  velocity  and  k  the  Boltzmann  constant. 

Different  linearization  or  expansion  techniques  are  appli¬ 
cable  in  cases  where  f  is  closely  approximated  by  this  Maxwellian  distri- 
ution,  the  best  known  being,  probably,  the  Chapman  Enskog  expansion  (ref.  2). 
This  technique  assumes  small  gradients  for  the  macroscopical  quantities 
(bulk  velocity,  density,  pressure  or  temperature)  and  lead  to  the  classical 
Navier-Stokes-Fourier  equations  as  a  first  approximation,  describing  the 
so-called  continuum  flow  regime. 


U 
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At  the  other  end  of  the  range,  the  free  molecule  regime, 
the  collision  term  becomes  negligible,  and  equation  (1.1)  reduces  to  a  near 
trivial  first  order  differential  equation.  The  interaction  between  mol¬ 
ecules  and  solid  boundaries  then  represents  the  central  problem  but  no 
major  mathematical  difficulties  are  encountered. 

All  other  regimes,  between  the  two  extremes  quoted,  require 
a  full  treatment  of  the  Boltzmann  equation.  Usually,  a  distinction  is  made 
between  several  such  regimes:  slip  flow,  transitional  regime  and  near  free 
molecule  flow  (ref.  3). 

One  must  also  be  aware  of  the  fact  that  equation  (1.1)  only 
accurately  describes  monatomic  gases,  and  should  only  be  considered  as  a 
physically  simplified  model  for  the  complicated  behavior  of  the  technically 
important  diatomic  gases.  In  many  attempts  at  solving  eq.  (1.1),  one  goes 
even  farther  than  this  and  assumes  a  particular  interaction  potential  for 
the  molecules  in  order  to  simplify  the  collision  term. 

These  particular  potentials  are  sometimes  very  far  from 
realistic.  The  most  spectacular  simplification,  as  already  indicated  by 
Maxwell,  occurs  when  a  repulsion  proportional  to  the  inverse  of  the  fifth 
power  of  intermolecular  distance  is  chosen  (ref.  2).  This  Maxwell  molecular 
model  is  not  realistic  either. 

However,  it  appears  that  no  physical  model  will  lead  to  a 
collision  term  simple  enough  to  allow  an  exact  solution  to  be  obtained  in 
the  general  case  or  in  non  equilibrium  cases  of  practical  interest. 

Many  important  mathematical  questions  appear  in  connection 
with  equation  1,  such  as  its  behavior  in  the  vicinity  of  free  molecule 
regime,  the  validity  of  expansion  procedures  and  of  several  iiiq>ortant 
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approximative  methods*  Such  methods  are  usually  introduced  without  attempt 
ing  rigorous  mathematical  justifications  (see  for  instance,  Mott  Smith’s 
treatment  of  the  shock  wave  in  ref,  4  and  Lee's  discontinuous  distribution 
function  technique  in  ref,  5), 

It  seem.s,  therefore,  that  a  mathematical  model  without 
immediate  physical  signif icance^  but  which  would  present  most  of  the  main 
features  of  eq.  1,1,  would  be  valuable  if  exact  solutions  could  be  obtained 
for  significant  cases. 

First  of  all,  one  should  know  what  are  the  most  significant 
features  of  equation  1.1,  from  the  mathematical  standpoint.  This  is  not 
easy  because  so  little  is  known  about  the  properties  of  the  solutions. 

The  mathematical  model  studied  in  this  report  presents  the 
following  basic  properties  in  common  with  eq.  1.1: 

a)  It  has  the  same  structure,  involving  space  and  time 
derivatives  in  the  first  member  and  a  "collision  term"  in  the  second  member 

b)  There  are  collisional  invariants  as  in  equation  1  (the 
number  of  molecules,  the  momentum  and  the  kinetic  energy),  so  that  a  set 
of  hydrodynamical  equations  may  be  derived  in  the  usual  way, 

c)  The  only  equilibrium  solutions  are  also  of  Maxwellian 
type,  as  given  by  eq.  1.2. 

d)  Though  no  property  analogous  to  Boltzmann’s  H  theorem 
(ref.  1)  could  be  obtained  for  the  model  equation,  it  appears  from  its 
analogy  to  the  Fokker  Planck  equation  (ref.  6)  that  the  collision  provides 
a  dissipative  mechanism. 
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e)  The  boundary  conditions  that  have  to  be  imposed  in  order 
to  define  a  solution  are  the  ones  that  would  be  physically  encountered  for 
an  actual  gas  kinetic  system. 

f)  It  can  be  shown  (see  paragraph  4)  that  the  model  equation 
exhibits  an  infinite  number  of  relaxation  times,  just  as  the  Boltzmann 
equation  (a  condition  that  is  not  satisfied,  for  instance,  by  the  simple 
Krook  model  of  ref,  7). 

These  appear  to  be  the  widest  possible  similarity  to  the 
Boltzmann  equation  one  can  hope  to  obtain  using  a  reasonably  simple  mathe¬ 
matical  model  (as  distinct  from  a  physical  model)  . 

One  of  the  difficulties  one  might  fear  because  of  the 
absence  of  direct  physical  meaning,  is  that  solutions  for  f  corresponding 
to  reasonable  boundary  or  initial  conditions  would  turn  out  to  be  meaning¬ 
less  because  f  would  take  negative  values  in  some  regions  of  the  field* 
or  exhibit  singularities.  The  same  difficulty,  of  course,  could  arise  in 
applying  approximate  methods  which  are  hot  mathematically  justified  to 
equation  1.1. 


However,  the  analogy  with  a  Fokker  Planck  equation  again 
enables  one  to  hope  that  negative  f  values  will  not  occur  in  solutions 
corresponding  to  physically  acceptable  boundary  and  initial  conditions 
(i .e .  conditions  defining  entirely  the  solution  and  such  that  the  initial 
and  boundary  values  of  f  are  everywhere  positive) . 

It  is  possible,  of  course,  that  the  positive  definiteness 

*  It  is  clear  that  such  difficulty  could  never  occur  in  a  small  disturbance 
treatment  where  f  *  f^  +  ^  f  and  f^  is  always  positive  and 
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of  f  could  be  proved  for  the  model  and  for  all  valid  boundary  and  initial 
conditions  but  such  a  propf  has  not  been  attempted. 

Attempts  will  be  made  to  obtain  exact  solutions  for  boundary 
and  initial  conditions  of  physical  interest  such  as  steady  heat  conduction, 
wave  propagation,  steady  shock  waves,  etc.  The  behaviour  close  to  the  free 
molecule  regime  and  the  validity  of  the  usual  approximation  methods  will 
then  be  studied. 


2 •  The  One  Dimensional  Model 

An  important  simplification  of  eq.  1.1  would  arise  if  a 
fully  one  dimensional  situation  could  be  considered.  This  would  lead  to 
the  physical  model  of  a  system  of  molecules  restricted  to  move  on  a  straight 
line,  say,  the  axis.  The  number  of  independent  variables  would  then  be 
reduced  from  seven  to  three  in  the  general  case  and  from  six  to  two  in 
stationary  situations. 

However,  such  a  model  is  useless  if  a  physical  collision 
process,  with  the  usual  conservation  of  momentum  and  energy,  is  assumed. 
Simple  exchange  of  velocity  would  then  occur  for  each  collision  (  ^ 
and  =  Si  )  SO  that  the  collision  term  would  be  identically  zero  (a 

one  dimensional  gas  is  therefore  always  in  the  free  molecule  regimel). 

The  relaxation  of  one  conservation  condition  could,  perhaps, 
save  the  one  dimensional  model.  However,  with  the  energy  condition  dropped, 
one  cannot  expect  a  Maxwell  type  equilibrium  distribution  to  exist,  while 
on  the  other  hand,  to  keep  the  conservation  of  energy  and  drop  the  conserva¬ 
tion  of  momentum  does  not  seem  to  lead  to  clear  physical  interpretations. 


?• 


It  thus  appears  that  no  simple  one  dimensional  physical  model 
can  be  imagined* 

The  mathematical  model  proposed  here  is  suggested  by  the 
study  of  the  Boltzmann  collision  terra  for  weak  collisions.  This  term,  which 
turns  out  to  be  quite  comparable  to  the  Fokker  Planck  interaction  term  used 
in  the  treatment  of  ionized  gases  (ref.  6),  is  mentioned  in  ref.  1  without 
details  or  reference  to  other  papers . 


This  term  is  therefore  calculated  in  appendix  I* 


Initially,  it  was  hoped  to  use  the  weak  interaction  term 
and  some  further  simplifications  (involving  quite  eccentric  interaction 
laws,  as  shown  in  appendix  I)  in  order  to  obtain  a  simple  physical  model. 
The  resulting  equation  is  still  very  complicated,  as  shown  below 


where  %  is  a  parameter  depending  on  the  law  of  interaction,(r|^^j(jthe  stress 
tensor  and  is  the  Laplace  operator  in  the  velocity  space.  ^ 

denotes  the  random  velocity  of  the  molecules: 


This  equation  is  again  a  non  linear  integro-diff erential 
equation  for  f,  but  the  integral  aspect  is  simplified  by  the  fact  that 
only  a  finite  number  of  '’moments”  appear  which  involve  integrals  over  f, 
namely  the  ten  moments 

ti^  Iff  f 


am 
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We  have  then  tried  to  obtain  a  one  dimensional  mathematical 
model  of  similar  structure  namely 


(2.2) 


where  A,  B,  C  must  be  functions  of  ^  and  simple  moments  such  as  n,  u  and 

The  only  significant  difference  with  (2.1)  appears  to  be 
the  acceptance  of  a  first  derivative  with  respect  to  the  velocity.  It  is 
not  clear  whether  this  has  a  particular  physical  or  mathematical  signifi¬ 
cance  . 


We  now  require  conditions  (b)  and  (c)  of  paragraph  1  to  be 
satisfied,  i.e.  that  eq.  2.2  would  satisfy  the  conservation  of  number  of 
particles,  momentum  and  energy  and  that  the  equilibrium  solutions  (i.e. 
solutions  of  the  equation  obtained  by  equating  the  collision  term  to  zero) 
would  lead  to  a  Maxwellian  distribution  (1.2).  As  shown  in  appendix  II, 

A,  B  and  C  are  entirely  defined  by  these  conditions  up  to  a  common  factor. 
Logically,  this  factor  could  be  a  function  of  n  and  (not  of  u,  of 
course,  because  of  the  principle  of  Galilean  relativity) . 

Physically,  one  iraist  expect  that  a  collision  term  ,  involv¬ 
ing  pairs  of  molecules,  should  be  proportional  to  the  square  of  the  density. 
The  factor  is  thus  written  as  Xfl/  ,  where  X  may  be  a  function  of  . 

For  simplicity  we  will  take  a  constant  value  for  this  parameter.  The 
resulting  equation  then  becomes 


(2.3) 
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with 
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"X  =  Const. 


^Oo 


+<\? 

rur^*/ 


(2.4) 


This  is  the  mathematical  model  we  intend  to  study.  The 
interaction  term  is  very  similar  to  the  Fokker  Planck  term,  except  that 
the  coefficients  of  the  derivatives  involve  integrals  of  the  distribution 
function  which  greatly  complicate  the  mathematical  nature  of  the  equation 
by  making  it  non  linear  and  integral. 


3 .  Basic  properties  of  the  model  equation  2.3 

Equation  2.3  has  been  obtained  under  the  conditions  that 
the  number  of  molecules,  momentum  and  energy  are  collisional  invariants. 
This  means  that  the  "collision  term"  has  no  global  effect  on  total  density, 
momentum  and  energy.  Integration  over  the  velocity  ^  ,  of  equation  2.3, 
when  multiplied  successively  by  ,  and  yield  the  usual 

macroscopic  equations  of  continuity,  momentum  and  energy. 
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where  the  density  p  ,  the  pressure  p  and  the  heat  flux  q 
are  defined  as  usual  in  kinetic  theories: 


p 


iOO 

^  ina= 


(3.4) 


p  -  ^  'wij  { 


(3.5) 


— 


(3.6) 


Equations  (3.1),  (3.2),  (3.3)  are  entirely  analogous  to  the 
system  of  equations  describing  the  one  dimensional  motion  of  a  hypothetical 
gas  with  a  ratio  of  specific  heats  ^  equal  to  three*  and  without  viscosity. 
It  is  clear  that  viscosity  cannot  occur  in  purely  one  dimensional  molecular 
motions . 


Temperature  may  be  defined,  of  course,  through  the  funda¬ 
mental  relation  of  statistical  mechanics: 


4  iT 


average  kinetic  energy  of 


4  7 


me 


so  that  the  usual 


random  motion  for  a  molecule 
equation  of  state  for  a  perfect 


gas  applies: 


The  only  deviation  from  ideal  fluid  flow  results,  in  the 
absence  of  viscosity,  from  the  presence  of  the  heat  flux  term  in  the 
energy  equation. 


*  This  is,  of  course,  the  value  to  be  expected,  from  statistical  mechanics, 
for  a  gas  with  molecules  having  only  one  single  degree  of  freedom  (i.e. 
translation  along  the  x-axis) . 
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The  macroscopic  equations  (3.1),  (3.2)  apd  (3.3)  are  studied 
in  more  details  in  appendix  III,  where  the  corresponding  Navier-Stokes- 
Fourier  equations  are  derived  in  the  usual  way,  as  a  first  approximation 
for  small  gradients.  Also,  the  existence  of  a  shock  wave  is  shown,  and 
its  fine  structure  analyzed  according  to  these  Navier-Stokes-Fourier  type 
equations , 


In  appendix  II,  the  condition  that  the  equilibrium  solutions 
be  of  Maxwellian  type  was  used  to  obtain  a  necessary  condition  for  the 
coefficients.  We  must  now  analyze  whether  this  condition  is  sufficient 
to  ensure  that  these  are  the  only  equilibrium  solutions. 


The  equilibrium  solutions  must  satisfy 


and  the  general  solution  of  this  ordinary  linear  differential  equation  for 
f  is  immediately  obtained  as 


C,  e 


c,  e 


f _ 

iF  I  0  Al' 


and  being  arbitrary  constants. 

The  second  term  is  of  course  a  particular  solution  of  the 
non  homogeneous  equation 

if  f  f  ^  Cj 

which  is  readily  seen  to  be  unacceptable.  Indeed,  f  must  tend  to  zero 
for  large  values  of  j'^jand  so  must  .  Therefore,  f  would  behave 

asymptotically  like  C'z  for  large  and  the  momentum  and  total 

energy  would  diverge . 
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We  may  thus  conclude  that  all  meaningful  equilibrium  solu¬ 
tions  of  eq.  (2.3),  yielding  finite  density,  momentum  and  energy  must  be 
of  Maxwellian  type. 


4 .  Mathematical  Nature  of  the  Model  Equation 

The  model  equation  (2.3)  unfortunately  exhibits  some  basic 
mathematical  complications  in  common  with  the  full  Boltzmann  equation  as 
well  as  with  some  other  models  extensively  used  (such  as  the  Krook  equation 
of  ref.  7),  namely  in  being  non  linear  and  integral. 

However,  one  can  imagine  a  procedure  for  solution  which 
formally  bypasses  these  difficulties,  because  they  both  arise  from  the 
appearance  of  the  parameters  yi  ^  M,  ^  C*'  in  the  coefficients  of  the  equa¬ 
tion. 


If  one  observes,  then,  that  equations  (3.1)  and  (3.2), 
resulting  from  (2.3),  provide  two  equations  for  these  parameters,  one  can 
visualize  the  following  formal  procedure  for  a  solution. 

Assume  a  given  set  of  initial  and  boundary  conditions,  suf¬ 
ficient  to  define  the  solution  f.  The  character  of  this  set  is  of  course 
defined  by  the  mathematical  nature  of  equation  (2.3),  considered  as  a 
partial  differential  equation  for  f,  whatever  the  nature  of  the  functions 
or,  ,  ^  being  essentially  positive.  The  linear  partial  dif¬ 
ferential  equation  in  f  can  then  be  solved  for  arbitrary  functions  C*' 

and  the  resulting  function  f  will  be  a  functional  of  H,  H, 

Tl  *  X  S  3 
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we  get  a  functionsl  equation  which  is  expected  to  completely  define  U;  6^ 
when  coupled  with  the  two  differential  equations  (3.1)  and  (3.2), 

One  could,  of  course,  also  devise  an  iterative  scheme, 
whereby  n,  for  instance,  is  chosen  arbitrarily,  or  better,  realistically, 
calculated  from  (3,1)  and  (3,2)  and  the  resulting  set 
substituted  in  (2,3)  and  the  resulting  linear  partial  differential  equation 
solved  for  f  taking  into  account  the  boundary  conditions.  00.^  are 

then  recalculated  from  this  function  and  the  process  repeated  indefinitely. 
In  order  to  justify  these  procedures,  one  must  of  course  analyze  the  nature 
of  the  functional  relationship  between  f  and  C,>  in  the  first  case, 

or  the  convergence  in  the  second  case.  This  has  not  yet  been  attempted* 

However,  if  one  assumes  at  least  one  of  the  precedures  to 
be  valid,  the  study  of  the  boundary  and  initial  conditions  pertinent  to 
the  problem  may  be  made  for  equation  (2.3)  considered  as  a  linear  partial 
differential  equation  for  f,  assuming  n ,  xt  ;  to  be  knoTO  functions  of 
X,  t  (with  the  restriction  that  riand  should  always  be  positive)  .  If 
these  functions  were  constants,  equation  (2o3)  would  reduce  to  the  standard 
Fokker-Planck  equation,  and  it  is  therefore  likely  that  the  initial  and 
boundary  conditions  suggested  by  the  physical  origin  of  the  equation  will 
apply.  These  are,  of  course,  the  same  as  would  be  considered  natural  for 
equation  (2.3). 


The  question  of  boundary  conditions  is  analyzed  in  more 
detail  ,  for  the  steady  state  case  in  the  next  paragraph. 


5 .  Mathematical  Nature  of  the  Model  Equation  for  Steady  State 

For  steady  states  (f  independent  of  t),  equation  (2,3) 

reduces  to 
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(5.1) 

Furthermore,  (3.1)  and  (3.2)  simplify  to 

2M^=0  err  -nU  -  A.  (5.2) 

30l 

^  Cn{ ^  enr  jrt  fu^+  C^)  (5.3) 

9CC 

Hence  mA  represents  the  mass  flow  and  mB  the  total  x-impulse  (momentum  flow 
plus  pressure) . 


Finally,  equation  (2.3)  becomes,  in  the  steady  state 


with  the  auxiliary  condition 


a.  /  fdl 


(5.5) 


The  constants  A  and  B  must  be  determined  after  solution  by 
solving  the  equations (for  some  particular  abscissa  x^) 5 

A  -  /  ,  3= 


This  is  a  set  of  two  equations  for  the  two  unknowns  A  and  B. 
If,  in  view  of  the  formal  procedures  described  in  paragraph  4,  one  may, 
for  a  moment,  consider  n  as  a  given  function  of  x,  the  equation  (5.4)  is 
seen  to  be  a  partial  differential  equation  of  parabolic  type.  However,  it 
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differs  from  the  simplest  and  classical  equation  of  this  type,  the  heat 
equation,  by  the  very  important  fact  that  the  sign  of  the  coefficient  of 

changes  in  the  field  (for  we  have  of  course  always  to  consider  the 
full  range  (^-oo  ^  j~oo)o£  values  of  ^  )*.  Gevrey  (ref.  8)  who  appears  to 
be  the  first  (and,  to  our  knowledge,  the  only  one)  to  have  dealt  with  such 
parabolic  equations  involving  a  coefficient  with  variable  sign,  calls  them 
singular  parabolic  equations.  He  also  shows,  in  his  fundamental  paper, 
how  the  usual  boundary  problem  pertinent  to  the  classical  heat  equation  must 
be  modified  in  the  case  of  singular  parabolic  equation. 


Before  quoting  his  results,. let  us  first  stress  that  the 
space  variable  x  in  h’quation  (5 .4)  plays  the  role  of  the  time  variable  in 
the  heat  equation  (because  it  appears  in  a  first  derivative  only),  while 
the  molecular  velocity  ^  plays  the  role  of  the  space  variable  in  the  heat 
equation  (because  it  appears  in  a  second  derivative) . 

The  simplest  boundary  value  problem  for  the  heat  equation 
(which  corresponds  to  the  simple  infinite  "wall"  problem)  consists  in 
giving  the  values  of  the  unknown  function  at  some  initial  time  and  also 
at  two  boundaries  (both  sides  of  the  "wall”)  at  all  subsequent  times.  In 
the  case  of  eq.  (5.4)  this  would  correspond  to  giving  f(0,^)^and  f(Cl,  ^  ) 
for  two  values  of  ^  .  Here,  however,  the  normal  boundaries  for  ^  being 

±  OO  ,  the  boundary  conditions  become  a  condition  about  the  asymptotic 
behaviour  of  f  for  large  absolute  values  of  ^  ,  as  already  mentlonned. 

An  important  remark  will  immediately  show  the  difficulties 
Connected  with  the  singular  parabolic  equation:  the  heat  equation  cannot, 
in  general,  be  solved  backwards  in  time. 

*  In  order  to  keep  the  usual  meaning  of  the  parameters,  ^  H,  C*-  should  of 
course  always  be  chosen  positive. 
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Indeed,  if  values  for  the  unknown  functions  are  prescribed 
at  some  final  time,  as  well  as  values  at  two  space  boundaries  at  earlier 
times,  then  the  solution  will  be  found  to  break  down  immediately,  unless 
very  special  final  values  are  chosen.  This  can  be  shown  easily  by  using 
a  Fourier  series  expansion  of  the  solution  in  terms  of  the  space  variable. 
The  coefficients  of.  the  expansion  are  found  to  be  exponentials  in  the  time 
variable,  increasing  for  decreasing  time  and  with  an  exponent  increasing 
very  rapidly  with  the  order  of  the  harmonics  (the  solution  will  therefore 
exist  in  a  finite  region  of  decreasing  times,  only  if  the  number  of  initial 
harmonic  components  is  finite  or  if  the  amplitude  of  the  harmonics  decreases 
sufficiently  rapidly  when  their  order  is  increased) .  This  phenomenon  is, 
of  course,  connected  with  the  extreme  "smoothing"  that  heat  conduction 
produces  on  temperature  distributions. 

This  remark  shows  then^that  a  classical  boundary  value 
problem,  with  one  Initial  and  two  boundary  conditions,  cannot  be  solved, 
in  general,  for  a  singular  parabolic  equation,  because  the  phenomenon  just 
mentioned  would  certainly  occur  in  the  region  where  the  coefficient  of  the 
time  derivative*  becomes  negative,  which  is  of  course  equivalent  to  a 
change  of  sign  for  the  time  variable  or  a  permutation  of  past  and  future. 

Gevrey  shows  in  ref .  8  that  the  standard  boundary  value 
problem  for  singular  parabolic  equations  must^as  expected^be  modified  as 
follows:- 


In  the  region  of  positive  coefficient,  initial  values  (for 
some  time  t^^,  say)  must  be  given,  as  well  as  a  boundary  condition  for  a 


*  When  referring  to  this  coefficient,  we  always  assme  of  course  that  the 
coefficient  of  the  second  derivative  is  equal  to  one  or,  at  least,  posi¬ 
tive  . 


17. 


value  of  the  space  variable  t  >  .  In  the  region  of  negative  coef¬ 

ficient,  final  values  (for  some  t^  ^  t^^,  say)  must  be  given,  as  well  as 
a  boundary  condition  for  a  value  of  the  space  variable  for  t  t^ . 

The  main  problem  arising  in  the  solution  is  then  seen  to  be 
a  matching  problem  along  the  "singular*®  line,  where  the  coefficient  vanishes 
Indeed,  assuming  arbitrary  values  for  the  unknown  function  on  the  singular 
line,  one  can  obtain  the  solution  in  the  region  of  positive  coefficient  just 
as  in  the  classical  heat  conduction  problem  (because  the  initial  values  are 
given,  as  well  as  two  boundary  values)  and  the  solution  in  the  region  of 
negative  coefficient  can  be  obtained  similarly.  However,  we  have  to  satisfy 
a  matching  condition  along  the  singular  line,  because  the  normal  derivative 
for  both  solutions  just  described  must  be  equal  at  all  points  of  this  line. 
This  condition  can  be  written  as  an  integral  equation  involving  the  unknown 
values  on  the  singular  line  which  are  then  obtained  by  solving  the  integral 
equation . 


In  the  particular  case  of  the  model  equation  (5,4),  the  sing¬ 
ular;  line  corresponds  to  ^  =  0  and  is  of  course  the  x-axis  in  a  X;  ^ 
diagram.  The  normal  boundary  conditions  are,  as  already  stressed,  replaced 
by  conditions  defining  the  asymptotic  behaviour  of  the  distribution  function 
This,  again,  may  be  considered  as  a  condition  for  ^  -  -HOO  and  another 
condition  at  ^  oo  ^  Xhe  boundary  conditions  for  the  Gevrey  problem, 
must  then  further  include  the  values  of  f  at  a  certain  abscissa  x^^,  for  the 
positive  values  of  the  velocity  ^  ,  as  well  as  the  values  of  f  at  another 

abscissa  for  the  negative  values  of  ^  . 

This  is  of  course  completely  in  agreement  with  what  would  be 

suggested  by  purely  physical  considerations.  Consider  a  finite ^linear 

volume"  of  one  dimensional  gas,  on  the  segment  say,  and  the  boundary 

conditions  at  the  two  solid  "walls"  at  x  ®  x,  and  x  *  x^ .  Interaction 

1  2 
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between  the  molecules  and  the  solid  walls  may  be  very  complicated,  but  will 
always  lead  to  an  expression  of  the  velocity  distribution  of  the  outgoing 
molecules  in  terms  of  the  distribution  of  the  incoming  molecules*.  There 
is,  however,  no  mechanism  described  in  the  kinetic  equation  or  in  the 

boundary  conditions  which  enables  the  walls  to  control  the  velocity  of 

a 

the  incoming  molecule,  for  this  would  involve  action  atYaistance.  Accord¬ 
ingly,  the  conditions  introduced  by  wall  effects  only  affect  outgoing 
molecules,  so  that  for  the  wall  on  the  left,  at  x  »  we  will  obtain 
information  about  the  distribution  function  of  the  molecules  with  ^  >  O 
and  at  the  wall  on  the  right,  at  x  =»  x^,  information  about  the 

distribution  function  of  the  molecules  with  .  This  is  similar  to 

a  Gevrey  problem.  Of  course,  the  situation  is  much  more  complicated  here, 
in  that  we  must,  in  general,  consider  relationships  between  the  values  of 
f  for  positive  and  negative  ^  at  the  walls,  and  only  in  the  case  where 
complete  accommodation  is  assumed  (i.e.  outgoing  molecules  have  no  memory 
of  their  condition  at  arrival)  do  we  get  the  simple  Gevrey  type  boundary 
conditions . 


We  may,  however,  conclude  that  the  model  is  very  satisfac¬ 
tory  from  the  point  of  view  of  the  nature  of  the  pertinent  boundary  value 
problems  for  the  steady  state,  which  are  in  complete  agreement  with  the 
physical  nature.  This  again  could  be  expected  from  the  analogy  with  the 
Fokker-Planck  equation. 


6*  Unsteady  Spatially  Uniform  Solutions  for  Equation  (2.3). 

It  is  interesting  and  comparatively  easy  to  examine  the 
solutions  of  (2.3)  depending  only  on  t  and  ^  but  not  on  x,  i.e.,  the 


*  For  instance,  through  an  integral  relationship  involving  a  stochastic 
kernel,  as  described  in  ref.  1. 
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solutions  uniform  in  space.  This  will  enable  us  to  show  that  the  model 
equation  exhibits  infinite  number  of  relaxation  times  and  modes.  The 
problem,  in  itself,  has,  of  course,  little  physical  significance,  for  the 
corresponding  initial  conditions  (i.e.,  a  non  Maxwellian  initial  distribu¬ 
tion  of  f,  uniform  in  space)  cannot  easily  be  created  and  does  not  occur 
in  practical  cases. 

The  analogous  problem  has  been  treated,  for  the  Boltzmann 
equation,  and  for  the  Maxwellian  molecule  model  by  Maxwell  (ref.  9)  and 
also  involves  an  infinite  number  of  relaxation  modes. 


In  the  spatially  uniform  case,  the  macroscopic  equations 
(3.1),  (3.2)  and  (3.3)  reduce  to 


-  c? 

JJc 


? 


and  ^ 


so  that  n,  u  and  C*-  are  constant  and  equal  to  their  initial  values,  and 
equation  (2,3)  thus  becomes  linear. 


We  first  introduce  the  non  dimensional  variables 


"\nt 


%  = 


fW 


<6.1) 


(reduced  time  and  random  molecular  velocity)  and  obtain  the  following  sim¬ 
plified  version  of  (2.3); 


2f 

9z 


(6,2) 


while  the  auxiliary  conditions,  defining  u  and  as  integrals  over  f,  may 
be  written 
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—  Qo 

T 

J  %  ^ 

(The  equation  defining  n  is  irrelevant  because  the  solution  f  will  always 
involve  an  arbitrary  constant  factor  which  can  be  adjusted  to  obtain  the 
correct  density) . 

There  are,  of  courses,  boundary  conditions  for  ^=±oo 
in  that  f  must  tend  to  zero  for  large  |  'E\  and  that  the  momentum  and 
energy  integrals  must  converge. 

The  resulting  mathenjiitical  problem  turns  out  to  be  entirely 
analogous  to  the  problem  of  computing  the  wave  functions  for  the  quantum 
mechanical  linear  oscillator  (ref .  10) . 


(6.3) 


(6.4) 


The  problem  is  treated  in  a  different  way  in  appendix  IV, 
using  a  Fourier  tran3forma.tion  with  respect  to  the  velocity  variable  ^  , 
because  this  technique  proves  to  be  useful  for  the  study  of  the  linearized 
model  equation  (see  also  paragraph  7s  eq.  7.15). 


The  general  solution  of  the  problem  is  given  by 

j .  e'^  / -e  a,  ■  j 


where  H  is  the  Hermit  e  polynomial  of  order  n  (ref,  11) 
n 


MiC?)  -  C"  -1 


etc . 
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The  auxiliary  conditions  must  now  be  taken  into  account. 
Multiplying  the  solution  obtained  by  and  (  )  respectively  and  in¬ 
tegrating  between  -oo  and  yield,  because  of  the  ortho¬ 

gonality  properties  of  the  Hermite  polynomials  (ref.  11): 


f  f  ‘  '/ijr  a,  6  ^ 

—  oO 

wpo  ~ 


The  auxiliary  conditions  therefore  imply 


The  other  coefficients,  etc.  are  well  defined 

by  the  Initial  condition 


t-o 

~?L 

because  of  the  fact  that  the  functions  C  ^  form  a  complete 

orthogonal  set. 


Hence,  it  is  clearly  seen  that  the  function  f  tends  to  a 
Maxwellian  distribution. 


i  - 


i'Ci 


for  t  ->■  oo  ,  and  this  tendency  is  characterized  by  an  infinite  set  of 
relaxation  modes,,  each  of  which  has  its  own  relaxation  time.  The  relaxation 
times  are,  in  decreasing  order 

" sTn,  ’  ‘tTn  ’  ‘  '  “<=■ 

and  correspond  to  faster  and  faster  relaxation.  The.  relaxation  times  are, 
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as  expected,  inversely  proportional  to  the  number  density. 

It  is  also  seeiig  from  the  form  of  the  solution  obtained  in 
appendix  IV,  that  f  remains  positive  definite  if  it  was  initially  positive 
definite . 


7  •  Remarks  on  the  Linear i  Blquation  (.2  *3) 

It  is  useful  to  attempt  the  solution  of  linearized  version 
of  complicated  equations  before,  attacking  the  fully  non  linear  case^s.  In 
the  case  of  the  model  equation  (2.3)5  i^^  case  of  the  Boltzmann 

equation  the  simplest  example  of  a  spatially  non  uniform  linear  problem 
occurs  when  considering  a  smils  steady  disturbance’' of  an  equilibriums 
Maxwell  distribution^  of  the  type 


may  be  a  spati.aliy  non  unifom  fuxiction  if  Tlc,  and  are:  considered 
to  be  functions  of  x  (loc.^il  fexweillan  distributiian.  See*,  for  Instauceg 

appendix  III)  . 

We  ass-amsy  here,  thsit  71^^  Uo  and  are  constants  so  that 
we  study  disturbances  from  an  absolute  evpilibrium  (uniforia  initial  condi¬ 
tions)  . 

In  this  c.'ase  we  writep  as  is  usual 

where  f*y  n®  and  S!xe  disturbancesp  which  must  be  considered  uniformly 

small  compared  to  the  corresponding  equilibrium  values  f  3  n  and 

c  o 

while  the  non  uniform  bulk  velocity  component  u®  is  limited  by  the  condition 
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u'<cfF 

a  condition  which  may  be  expressed  by  saying  that  the  Mach  number  based  on 
the  velocity  disturbance  must  be  very  small. 


The  linearized  version  of  equation  (2.3)  is  obtained  by 
neglecting  the  terms  of  order  higher  than  one  in  the  disturbances s 


Equations  (3 .1)  and  (3 .2)  become 


n 


W  - 55 


(7.2) 


(7.3) 


If  we  consider  the  steady  state  case,  the  time  derivative 
drops  and  the  equation  (3.1)  and  (3.2)  give 

(jliA)' ^  M  (7.4) 


(7,.5) 


It  is  clear,  however,  that  by  a  suitable  choice  of  the 
reference  values  Cp*'  ,  one  can  make  A'  and  B'  vanish,  which 

yields 


~  -f-  yi ' [Av  ~  o 

(uu)'Up 
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so  that,  with  this  choice  of  reference  values 


ft' 

rio 


(7.6) 


(fic^y^-ncK'U^^  ut  n' 


(7.7) 


It  is  useful  to  introduce  the  non  dimensional  variables 

S  not 


'  ,nr 


(7.8) 


as  well  as  the  parameter  N*,  related  to  n*  by 

400 


!\J'-  ^ 


(7.9) 


/ 

-oa 

The  reduced  equation  then  becomes 

Ft  *  (7.10) 


where 


is  a  non  dimensional  parameter,  given  by 


t^o  ^  ^ 

^  iW 


(7.11) 


which  is,  of  course,  closely  related  to  the  Mach  number  ^M= 

In  problems,  such  as  heat  transfer,  where  no  bulk  velocities 
exist,  we  must,  of  course,  take  so  that  the  equation  is  slightly 

simplified  to 


e  ^ 
fix 


(7.12) 
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The  Fourier  integral  transformation  with  respect  to  the 
(reduced)  velocity  component  ^  ,  as  used  in  appendix  IV,  also  provides 

a  simplification  in  the  actual  case.  Let  (p  be  the  Fourier  transform  of 
f which  exists  in  view  of  the  fact  that  f  must  tend  to  zero  for  large 
values  of  \%\  in  the  case  of  physically  valid  solutions.  We  have 

+  90  ■  ■  4^  ^ 


We  then  transform  the  equation,  taking  into  account  the  fundamental  pro¬ 
perties  of  the  Fourier  transformation  (ref*  12) 

‘-r  s/  ^  (7.U) 


and  obtain 


(7.15) 


Using  the  obvious  formula 


N'=  /  /'<*?  =^(p(r),X,o) 

-pc 


(7.16) 


and  the  fact  that  the  function  0.  ^  is  invariant  under  the  Fourier 

transformation  (ref.  12). 


In  the  case  of  a  steady  state,  it  is  seen  that  the  singular 
parabolic  equation  is  reduced  to  a  hyperbolic  equation. 


The  simplest  physically  significant  problems  that  arise  in 
connection  with  these  equations  are  the  steady  state  heat  conduction  pro¬ 
blem  and  the  wave  problem. 
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In  the  steady  state  heat  conduction,  one  assumes  tip  -=  O 
(and  therefore  =  O  ),  so  that  the  following  direct  and  trans¬ 

formed  equations  have  to  be  used 


hJiffdZ  0.17) 


(7,18) 


If  the  heat  conduction  is  treated  for  a  finite  amount  of  gas 
between  two  walls,  at  ^  and  1  ^  correspond^ 

ing  boundary  conditions  would  be  of  the  type 


£or 

(7.19) 

=  for  '^  <  CP 


Usually^  one  would  assume  that  the  molecule  ref lected  from 
the  wall  are  completely  accomodatedj,  thermally,  with  the  wall.  If  the 
reference  values  ^  are  chosen  to  correspond  to,  say,  the  conditions  at 
the  wall  at  ^  ,  we  would  then  have  the  simplest  set  of  possible  boun¬ 

dary  conditions 


for  >  <9 

for  <  O 


(7.20) 


The  value  for  is  arrived  at  by  writing  the  Maxwellian 

distribution  ^  corresponding  to  the  parameters  C/"- 

and  computing  the  first  order  term  of  f^  -  f^,  taking  into  account  the  (in¬ 
tegrated)  macroscopic  equation 


01 2, 

The  treatment  of  the  equations  for  linearized  heat  conduc¬ 
tion  still  appears  to  be  very  involved.  Results  obtained  in  this  direction 
will  be  published  shortly. 


In  the  wave  problem,  one  investigates  unsteady  sinusoidal 
solutions  of  the  form 


CvV-(- 


(7.21? 


where  (j  is  an  amplitude  function.  This  is  expected  to  lead  to  an 

eigenvalue  problem,  yielding,  for  each  value  of  the  frequency  luf  ^  to  an 
infinite  set  of  propagation  modes .  We  have,  of  course 


(7.22) 

where  a  is  the  speed  of  sound  for  a  particular  mode  and  (p<  the  corresponding 
attenuation . 


The  equation  for  becomes  (assuming  the  medium  at  rest, 

and  hence  -  O  ) 


fee  0-X  ) 


fm 


O  (7.23) 


and  the  transformed  equation  simplified  to  a  first  order  equation 


Co' 


where 


is  the  Fourier  transform  of  g. 


(7.24) 


The  reduction  in  order  of  the  equation  from  second  to  first 
brought  about  by  the  transformation  is  of  course  very  useful,  for  the 
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solution  of  the  latter  equation  can  be  obtained  at  sight.  It  is  therefore 
interesting  to  note  that  this  reduction  corresponds here,  to  the  rejection 
of  (physically)  meaningless  solution  types.  We  already  noted,  in  paragraph 
3»  that  there  are  two  types  of  equilibrium  solutions,  one  of  them  being 
rejected  because  it  leads  to  infinite  total  momentum  and  energy.  It  is 
easily  seen  that  this  is  just  the  type  of  function  for  which  the  Fourier 
transformation  of  the  equation  is  notallowed,  because  terms  such  as  *2^ /• 
lead  to  diverging  transforms 

The  Fourier  transformation  thus  automatically  eliminates 
the  unacceptable  solutions  and  a  corresponding  simplification  of  the  trans¬ 
formed  equation  had  to  be  expected. 

The  simplifications  are  not  so  obvious  in  the  general  case. 
The  complete  treatment  of  the  wave  problem  will  be  presented  in  a  companion 
report . 


8 .  Preliminary  Remarks  on  the  Near  Free  Molecule  Regime 

One  of  the  most  interesting  and  difficult  problems  arising 
in  the  application  of  kinetic  theory  to  aerodynamic  problems  is  the  study 
of  the  near  free  molecule  regime^  or  in  other  words,  the  way  in  which  so¬ 
lutions  of  the  Boltzmann  equation  behave  for  very  large  Knudaen  number. 

The  use  of  the  model  to  clarify  the  situation  is  therefore 

indicated. 


Actually,  we  cannot  strictly  define  a  Kmidsen  number  for  the 
model,  because  the  collision  term  cannot  be  interpreted  as  representing  two 
body  collisions  and  the  notion  of  mean  free  path  therefore  disappears. 
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However,  it  is  clear  that  the  free  molecule  regime 
corresponds  (for  given  boundary  conditions)  to  a  vanishing  We  must 
therefore  study  the  behaviour  of  solutions  corresponding  to  given 
boundary  conditions  for  X  tending  to  zero. 


The  following  remarks  concern  the  model  equation  for 
steady  state  i.e. 


which  reduces,  for  A  =  0  to 

7'X. 

It  is  clear  that  the  type  of  the  equation  entirely  changes  for 
so  that  one  has  a  singular  perturbation  problem  for  small  A  . 


(8.1) 

(8.2) 

0, 


Let  us  consider  the  simplest  case  of  heat  conduction 
between  two  walls  at  x  =  ^1^  (u  =  0) .  The  physical  boundary 

conditions  are  known  to  be  of  the  type  (Cf  7.19) 


for  •§  <  o 


and  the  corresponding  free  molecule  solution,  (for  A  *  0) ,  is  of 


course 


%>o 

I  %^o 

In  general,  we  therefore  expect  singularities  (discontinuities  of  ^ 
its  derivatives)  at  ^  =0.  The  approximation  for  the  free  molecule 
solution  must  therefore  break  down  in  the  vicinity  of  ^  =  0  in  this 
simplest  case.  A  small  region  will  exist  there,  say  from^  -  -  to 
^  +  t,  where  the  interaction  term  still  has  to  be  taken  into 


(8.4) 


or 


30, 


account.  The  situation  here  is  quite  similar  to  the  one  leading 
to  the  occurence  of  boundary  layers  in  fluid  flows  with  small 
viscosity,  where  viscosity  terms  still  have  to  be  taken  into  account 
in  regions  of  small  extension  %lose  to  the  walls  (as  well  as  in  the 
wakes) , 


One  can  readily  evaluate  the  importance  of  the  region 
where  the  interaction  term  is  important^  using  an  order  of  magnitude 
argument  similar  to  the  one  used  to  compute  boundary  layer  thickness 
in  terms  of  Reynolds  number  in  viscous  flow  theory. 


take  ^  ~  ^2 

to  be  of  order 
for  L  :  /—  = 


Let  L  be  a  length  scale  for  the  problem  at  hand.  We  would 
in  the  heat  conduction  problem.  We  then  assume 

(this  in  factj  leads  to  more  rigorous  definitions 
1  1 


max 


9  for  instance).  Let  L  be  the  width  of  the 


region  around  the  line  ^  =  0  in  which  the  interaction  term  is  important. 

The  definition  of  ^  also  has  to  be  made  more  precise^  by  stating  that 
9f  f  ^ 

is  to  be  of  order  ^  and  of  order  in  the  said  region. 

The  first  member  of  equ.  (8.1)  is  seen  to  be  of  order  ^  L.  and  the 


second  member  involves  terms  of  the  order  of  and  i 

the  region  between  ^  .  If  £  is  small  compared  to  the  velocity 

scale  for  the  problemg  the  first  term  of  the  second  member  is 

seen  to  be  dominant  and  an  order  of  magnitude  comparison  then 
follows  : 


n 


and,  hence 


so  that  £ 


—  ro  ^  n  c*'  ^ 

L 

£  oJ>  n  C’' 

is  small  of  order  ~  when  \  tends  to  zero.  Furthermore, 
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equ.  (8.1)  reduces,  in  first  approximation,  to 


3^^ 


(8.5) 


in  the  singular  region.  This  equation  is  analogue  to  the 
boundary  layer  equation  in  viscous  flow  theory.  It  is  seen  to  be 
equivalent  to  the  typical  singular  parabolic  equation  treated  by 
Gevrey  in  ref. 8,  because  reduces  to  a  constant.  It  is,  of  course, 

useful  to  introduce  suitable  non  dimensional  variables 

3:  =  ^  ^  -J _  X  -  ^ 

and  (8.5)  then  reduces  to  the  Gevrey  canonical  form 

—  9i- 

9X 


The  boundary  conditions  (8.. 3)  in  the  singular  region, 
become  for  \  tending  to  zero 


iCo,  ftrr  ^ 

j-C'i, ^  firr  <o 


(taking  L  ~  ^2  ~  1^  origin  of  the  x  axis  being  displaced  to 

the  point  x  =  x^^) . 


This  is  valid  for  (o)  If 

but  ^  ^CO)  ^  must  write,  in  the  singular  region 

/  - 


obviously  satisfies  the  same  equation  (8.6)  but  the  pertinent 
boundary  conditions  are  now 


II 

O 

00 

I 'I 

-frr  \—>C? 
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The  contribution  to  from  the  interaction  term 

is  therefore  of  order  /  in  the  singular  region  of  thickness^ ^ 

and  of  order  (r  elsewhere  (i  e  for  (C/  ^  £).  If  we  then 

1^1^ 

compute  the  influence  of  the  interaction  term  on  the  deviation  from 
free  molecule  values^  for  physically  significant  macroscopic  quantitieSj 

— — r  • 

such  as  cj  etc,  3  the  contribution  of  the  ^’boundary  layer®*  will 

be  of  the  order  )f/3 

and  the  contribution  of  the  other  regions  of 
order  A,  The  first  contribution  is  thus  seen  to  be  the  main  term 
in  deviation  from  free  molecule  values.  This  deviation  is  then 


proportional  to  j(/3 

for  0. 


The  details  of  the  calculations^  based  on  the 
solution  of  equ. (8«6)  with  the  boundary  layer  conditions  (8.7)  or 
(8,8)  will  be  presented  in  a  forthcoming  report. 
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APPENDIX  I  “  The  Boltzmann  Equation  for  Weak  Interactions 


We  write,  using  the  notation  defined  in  paragraph  1,  the 
momentum  and  energy  relations  for  a  collision,  between  molecules 
of  equal  mass  : 

•ICU  (I.l) 


Introducing  the  relative  velocities  S  and  <5*  with  respect  to  the 
center  of  mass  of  the  ro.olecule  pair  before  and  after  collision 


I*  =  2:^ 
2 


(1.2) 


we  get,  from  the  momentum  equation  : 

and  from  the  energy  equation 

I?'!  =151 

If  we  put 


0-5*4 


(1.3) 


(1.4) 


(variation  of  relative  velocity  caused  by  the  collision) ,  we  may 
use  the  following  series  expansion,  which  will  be  useful  to  analyse 
weak  interactions^  corresponding  to  small  values  of  |  ^  1 

and  a  similar  expression  for  ^(  *)  ^replacing  ^  by  ^  and  ^by  -(j) 


(1)  The  Einstein  summation  convention  is  used  throughout  this  appendix  : 
If  an  index  appears  once  in  a  term,  one  should  take  the  sum  of  the  terms 
corresponding  to  the  values  1,2 ,3  of  the  index,  unless  stated  otherwise  : 


therefore  means 
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After  some  algebra,  the  followixig  expansion  is  obtained 
for  the  main  factor  in  the  Boltzmann  collision  term 

'I  A  £> 'cffix)  ill 


.1 

G 


(1.5) 


Fig.  1  indicates  the  geometry  of  a  collision,  as  observed  from  the 
center  of  massjof  the  system  formed  by  the  two  colliding  molecules 
(thi^  center  of  mass  of  course  moves  with  uniform,  velocity  during  the 
\\7hole  collision  process,  because  of  the  principle  of  equality  of 
action  and  reaction).  Q  is  the  scattering  angle  and  we  have,  clearly 


(f)=  S^-6  =  SCccse 

where  ^  is  a  vector  perpendicular  to  £,  If  we  assume  weak 
interactions,  0  will  be  very  small  and  we  use  the  expansions 


(1.6) 


'tl-lSife-f,...) 


(1.7) 


Now,  the  collision  term  obviously  involves  an  integration 
with  respect  to  £,  the  angle  defining  the  position  of  the  plane  TT 
of  the  relative  motion  (relative  to  a  reference  plane) .  We  have, 
by  symmetry  : 


j  ^  cie  -Q 


which  leads  to  the  important  conclusion  that  no  term  of  first  drder 
in  Q  is  going  to  appear  in  the  collision  integral* 


In  order  to  compute  further  terms  of  this  integral,  we 
need  to  integrate  products  of  components  of  ^  with  respect  to 
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It  is  clear  that  many 


A  mathema^tical  trick  greatly  simplifies  these  laborious 
computations.  Let  CL  be  an  arbitrary  vector,  OC  the  angle  between 
the  vectors  and  0  and  S  the  angle  between*^  and  the  projection 
of  (X  on  a  plane  perpendicular  to  <5*  (see  fig  1) 

We  may  choose  the  arbitrary  origin  for  the  angle  to  coincide 
with  this  projection  so  that  S  =  ^  and 

j^l  9vw.o<  1'^ I  CoS  8  ^  1^1  yuA.oi'l  d  j  ^ecosC  —  I  J  |  ^O-CasS 


Th(g  method  suggested  for  the  simplification  of  the 
integrals  consists  in  using  the  following  identity,  for  a  product 
of  n  components  of  (some  of  which  may  be  identical) 


dajda^- 


a:c) 


From  this,  one  immediately  concludes  that 

The  integral  can  be  transformed  into  a  tractable  form  as  follows 


J 


(1) 

A)  .  ,  .11 

/.I 


I  yi  Z  (jxl) 


using  the  identity 


(a.Aa)‘^=  rV'-ai)'-. 


The  integral  can  actually  be  expressed  in  terms  of  Bessel 
functions  but  this  does  not  seem  useful  for  the  present  purpose. 
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Computation  of  the  integrals  over  of  and 


then  becomes  immediate 
,  ^ 


(1.8) 


il  =  o 


(  denotes  the  Kronecker  symbols  5^^'^  =  0  for  i  »  1  ior  ^  =J  ) 


i 


After  integration  with  respect  to  we  therefore  get  the  following 
expansion  in  powers  of  0  s 


A  ^ 

The  remaining  terms  are  of  order  y  or  higher 


(1^9) 


The  first  approximation  for  weak,  interaction  therefore 
would  lead  to  a  collision  integral  of  the  form 


(I. 10) 


(the  next  term  being  of  order  4) 

The  dynamics  of  a  collision  must  now  be  analyzed  in  more 
detail  in  order  to  obtain  the  relationship  connecting  the  scattering 
angle  ^with  the  collision  parameter  and  the  relative  initial 
velocity  |  ^  I  =  2^  .  In  the  spirit  of  the  weak  interaction 

assum.ptiong  one  should  again  try  to  obtain  an  expansion  in  powers  of  a 
parameter  indicating  the  order  of  magnitude  of  the  force  of  interaction 
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The  first  term  in  such  an  expansion  is  very  simple  for  it  can  be 
computed  assuming  that  the  trajectories  are  undisturbed.  The  final 
lateral  momentum  in  the  relative  motion^iTi^,  induced  by  the  central 
interaction  force  F  is  simply  obtained  as  the  total  impulse  produced 
by  this  force  in  the  direction  perpendicular  to  the  undisturbed 
trajectory;  i.e.  _ 

Now, 

I'd  =  I ^1  5i'n,0 

So  that,  in  the  first  approximation 

6  ^ 

This  leads  to  a  singular  factor  ^  integral  (1. 10)  which  is, 

however,  immediately  seen  to  converge  (because  of  the  first  factor  which 
involves  terms  of  order  — “^'1  and 
If  F  is  a  power  law  interaction 

Ifl-Kr'^  (■s>o) 

then  <\»  -sj„ 

/  F  4^  =  K 


(p=r/t) 


so  that  0  is  proportional  to  h'  ^  and  9^i}  proportional  to  ^  ^  . 

Therefore,  the  integral  fOhdb  would  converge  for  S^Z  and  diverge  for 

0 

^^2.  ^  the  second  case  certainly  showing  a  predominance  of  weak  inter¬ 
actions.  The  divergence  of  the  integral  has  to  be  remedied  by  considering 
a  finite  cut  off  distance,  as  is  done  in  the  theory  of  ionized  gases. 


We  may,  of  course,  consider  a  force  law  which  would 
nowhere  lead  to  strong  interactions,  for  instance  by  considering 
that  F  is  bounded  forf-^0. 
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The  collision  integral  (1,10)  is  still  very  complicated,  however, 
because  of  the  coupling  factor 

il-fl  fe^hih 

o  ^  -^1-3 

which  is  proportional  to  |  ^  .Maxwell  has  shown  that  this 

coupling  factor  disappears^  in  the  full  Boltzmann  collision  term, 
when  the  interaction  force  is  proportional  to  ^  (Maxwell" s 
molecular  model). 

It  would,  of  course,  be  meaningless  to  consi4er  weak  interactions 
for  such  a  power  law  in  view  of  the  obvious  predominance  of  the 
strong  interactions  in  the  collision  term. 


However,  one  could  imagine  a  rather  artificial  interaction 
model  which  would  induce  in  the  collision  term  (1,10),  the  same 
simplification  as  the  Maxwell  interaction  induces  in  the  Boltzmann 
collision  integral.  In  order  to  do  this,  one  must  assume  the 
interaction  force  to  depend  on  the  relative  velocity  of  the 
colliding  molecules;  i.e.  r  must  be  proportional  to  the  cube  of 
this  relative  velocity.  If  we  then  put 

^  i If-'C'I  X 

X  being  now  a  constants  integration,  with  respect  to  ^  “  becomes 
feasible  in  I.IO.  The  derivatives  of  ^ “  may  be  eliminated  by  partial 
integration  and  the  terms  at  the  limit  will  vanish  because  of 
the  required  asymptotic  behaviour  of ^ 

The  following  integrals  are  then  seen  to  appear  as 
of  ^  and  its  various  derivatives  in  1.10  i 

Iff  f' ^n. 

Iff  _ 

ill  /'  ‘‘■ti't  sy  ^  1' = ^ 


coefficients 
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Tlr  being  the  number  density,  lA.  the  bulk  (or  macroscopic) 
velocity,  XVj  the  stress  tensor  and  C  the  random  molecular  velocity. 

The  collision  integral  .for  this  model  reduces  to 

Xn  I  -h^J- }  (I. 

where 


is  the  Laplace  operator  in  velocity  space. 

The  collision  integral  therefore  appears  as  a  linear 

second  order  differential  form  in  ,  the  coefficients  of  which  are, 

in  fact,  integrals  over  j'  .  This  collision  term  is  accordingly 

still  non  linear  ^nd  integro  differential  with  respect  to  . 

However,  the  non  linearity  and  integral  dependance  here  appear 

through  a  limited  number  of  parameters,  namely  the  10  moments  iTiy 

and  which  might  produce  considerable  simplification  in  the 

J 

treatment  of  this  model. 


2. 1 


11) 
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APPENDIX  II  -  Derivation  of  the  Interaction  Term  Satisfying  the 


Condition  of  Paragraph  1. 

The  interaction  term  is  of  the  form 

where  Aj  Bj  C  should  be  functions  of  ^  and  simple  moments  of  the 
distribution^  such  as  n,  a-rid 

The  conservation  of  total  mass_, momentum  and  energy  in 
presence  of  an  interaction  of  this  type  requires  the  following 
integral  conditions  to  be  satisfied 

+  =o 

The  derivatives  of  j-  may  be  eliminated  by  partial  integration 
(assuming  all  integrals  to  exist)  c,  and  the  following  conditions 
for  A,  Bj,  C  are  then  obtained  i 

§CJ  -O 

II-2 

It  is  immediately  clear  that  A^  B^  C  must  depend  on  some 
moments  of  ^'^for  if  Ap  Bp  C  were  completely  independant  of  -f  ^ 
and  II. 2  was  to  be  satisfied  for  all  -j’  pone  would  be  led  to  a 
trivial  solution  A  =  B  =  C  =  0 

The  assumption  that  Ay  Bp  C  are  to  depend  only  on  riylTlUand  p  is 
equivalent' to  satis.fying  condition  II. 2  for  satisfying  the 
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auxiliary  conditions 

u-  Jfd? 

UK  =  j  f  ^  II. 3 

jD  ^  m 

yp  being  considered  as  given,  fixed  quantities. 


It  is  clear  that  these  three  conditions  reduce  to  two, 
because  the  equatbns  II.2  are  homogeneous  and  one  of  the  conditions 
II. 3  may  always  be  satisfied  by  adjusting  an  arbitrary  factor  ^  • 
We  must  then  reduce  II. 3  to  a  set  of  two  homogeneous  conditions, 
for  instance 


=o 


(  being  the  mean  square  of  the  random  velocity ^-6i) .  II. 2  must 
now  be  a  consequence  of  II. 4  for  all  satisfying  only  II. 4  .  We 
thus  have 


'i|  +  ^  =  «  (l-U.)  +<3(Xf5-<-c‘] 
fixrA_di!6  +  =  Yd-vi)  + 

df-  37^  " 

where  ^  are  independant  of 


It  is  useful  to  introduce  the  reduced  random  velocity 
^  "5  -U.  C 

{W  "(W 

and  to  combine  equations  II.  5  linearly  among  themselves  to  get  the 
reduced  system 
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4fe)  _  ^  C*  =  f}*?  +-  |S '*(?*- -()=  M 

4;^")_«l(B»r)  ,  r^c*  _  f  c  +  /'Y'?:'-1)  -  ^ 

((C‘  ‘d^ 

with 

Combining  the  equations  II.  6  linearly  among  themselves  (with 
respective  coefficients  ^  -l^and  1,  we  get 

2  A  -  N  „.7 


Multiplying  the  first  equation  by  ^  and  substracting  the  second  yields 


B*=  4L^-M 


II. 8 


and  finally 


II. 9 


We  must  now  return  to  the  third  condition  set  in  paragraph 
1,  namely,  that  the  Maxwellian  distribution  , 

j-  =  c<ov5t.e  =  c^h  e 

be  a  solution  corresponding  to  equilibrium,  or,  in  other  words,  that 
substitution  of  this  function  in  the  interaction  term  will  yield  zero. 
In  terms  of  the  unknown  coefficients  A,  B^,  C^,  this  reduces  to  the 
condition 


Ar-c'-i)  -3%  4  c*  =o 


II.  10 
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TSfow,  according  to  (II. 6)  (II. 7)  (II. 8)  and  (II. 9),  A,  B*,  and  are 
seen  to  be  polynomials,  the  degrees  of  which  are  respectively 
4,  3  and  2  at  most  (Uyt^/N  being  quadratic  in^  ).  According  to 
(II.  10),  however,  A  can  only  be  of  degree  2  at  most,  so  that  we 
must  have  (see  II. 7) 

L^O  or  =  0 

But  the  Last  result  shows  B*  to  be  a  linear  function  of  "2^ 

(see  II. 8)  and  we  conclude  then,  from  II.  10  that  A  must  be  a  constant 

A  =irN-2^M)=  coast. 

Introducing  the  explicit  expressions  for  N  and  H  then  yields 
and  hence 

We  have  then, finally 

A  =  const. 

B*  =  A 
C*  =  A 

These  values  obviously  satisfy  the  system  of  equations (II. 6)  and 
(11,10).  Reverting  to  the  original  variables  and  coefficients  we 
get,  for  the  interaction  term 


A  is  a  function  of  C*'  and  IX  which  is  taken  to  be  ^ 
in  paragraph  2. 
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APPENDIX  III  -  Hydrodynamic  Equations  Corresponding  to  the  Model 


The  following  macroscopic  equations  are  immediately 
obtained  when  the  model  equation  (2.3)  is  multiplied  successively 
by  and  ^i^d  integrated*  over  the  velocity  component^ 

from  -  00  to  +  CO  taking  into  account  the  conservation  propertj.es 
for  the  interaction  term  t 


Continuity  : 


^  ,  9 

2t  d%- 


Momentum  : 


Energy 


III.  1 


III. 2 

III.  3 


with  the  heat  flux 

—  oo 

If  we  consider  a  situation  for  which  the  deviation  from 
equilibrium  is  negligible^,  the  heat  flux  disappears  (for^  is  then 
a  Maxwellian  distribution)  and  we  obtain  the  equations  of  one 
dimensional  motion  for  ideal  gas  with  y  =  3s,  as  expected  from 
statistical  mechanics . 


As  already  stressed^  the  only  effect  of  a  deviation  from 
equilibrium  is  the  appearance  of  a  heat  flux<,  because  no  viscosity 
can  appear  in  a  hypothetical  gas  with  purely  one  dimensional  molecular 
motion. 

In  order  to  derive  from  (2.3)  a  complete  system  describing 
the  hydrodynamics  of  the  one  dimensional  gaSj,  similar  to  the  classical 
Navier  -  Stokes  -  Fourier  eq^uation%  a  technique  similar  to  the  various 
expansion  procedures  applied  to  the  Boltzmann  equation  will  be  used 
(Cf  Ref. 2) 


45, 


Obviously  we  require  a  formula  for  q  in  terms  of  the 
gradients  of  the  other  physical  parameters  (n,  u,  p) 

# 

To  obtain  an  equation  forq,  we  transform  equ. (2.3)  using 
the  variable  C  rather  than  This  yields 


^ 


(III. 5) 


1  3 

Multiplication  of  both  members  by  j  me  and  integration  with  respect 
to  C  produces  the  required  equation  of  transfer  for  q  .  Partial 
integrations  must  be  carried  out  on  the  interaction  term  in  order  to 
eliminate  the  derivative'5^'/^ with  respect  to  C.  The  result  is  : 


If  f 


We  now  look  for  a  first  approximation  to  this  relationship  under  the 
assumption  that  the  space  and  time  derivatives  of  are  small  and 
that  ^  is  well  approximated  by  a  local  equilibrium  distribution 
It  is  then  immediately  seen  that  ^  is  a  first  order  term  in  the 
derivatives,  (so  that  the  derivatives  of  and  ^  neglected) 

and  that 


inc^fdc 


ro 


neglecting  higher  order  terms. 
But 
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The  first  approximation  to  the  equation  of  transfer  for 
is,  accordingly,  given  by 

Q 


=  -5^'  tia 

I "  "  ^  j  1  i  f>  9xJ  ^ 


SO  that 


A  I  mc^  9c 


_  _  ^  m,  -r-^T 

I  'A  9%  z'X  d'%- 


(III. 6) 


(The  temperature  T  being  defined  by  Equ.  3.7  ) 

Formula  (III, 6)  is  equivalent  to  th^  Newton-Fourier  law 
of  heat  conduction,  the  heat  transfer  coefficient  being  given  for 
the  one  dimensional  model  gas,  by 

L  . 

1  X 

It  is  seen  to  be  proportional  to  the  temperature  and  in4ependent 
of  pressure,  just  as  in  the  case  of  the  Maxwell  molecular  model 
with  inverse  fifth  power  interaction. 


Equation  (2.3)  is  thus  seen  to  provide  a  dissipative 
mechanism  for  least  if  the  disturbance,  from  equilibrium 

is  small. 

Let  us  now  consider  some  simple  steady  state  problems 
for  the  model  equation 

In  the  steady  state,  equations  III. I,  2  and  3  reduce  to 


9  % 


p  u.  I  fA" + 1  £]  +  2£ir  =  -2i 

'  ^  2-  p/  92  9Z. 
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These  may  be  immediately  integrated  as  follows  (see  also  equ.(5.2)  (5.3) 

p  it  =  =  m  a 

-t-  p  =  ^  =  11x3 

being  constants) 


(III. 7) 


The  simple  heat  conduction  problem,  with  *1^  =  0  ,  is 
immediately  solved  as  follows,  in  the  first  approximation 

A  =  'Z  =  CifyiS-t 

q  ^ ylXs:  ^-{j-  con%C  ■ 

I  '  9% 

and  hence 


(III. 8) 


where  and  are  the  temperatures  at  the  "walls",  located  at  the 
abscissae  and 

The  ideal  fluid  approximation  (corresponding  to  ^  =  0) 
leads  to  the  classical  shock  condition,  but  for  the  case  ^  “  3 

p  btj  ^  =  7Kd 

(5u,^  +  p, 

If  we  introduce  the  Mach  number  M,  and  the  speed  of  sound  a  : 


(III. 9) 


M=  CL  =  u/y/^p/p 

the  following  relations  are  easily  derived  from  the  shock  conditions  : 

-:5  M|  — -j  =0  1 

-  i  fu, -f.  J  (III- 10) 
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The  detailed  structure  for  a  weak  shock  may  be  obtained 
approximately  by  using  the  first  approximation  described  by  the 
expression  (III06)  for  (j  and  substituting  this  in  equation  (III«7)o 
We  have 


-RT-  p/p  =  ^ 

Q^  —  fTnMX)  yc  —  (s/z)^‘U  -h 

'  dx. 


An  ordinary  differential  equation  for  li  is  readily  obtained  by 
eliminating  RT 


dn  ^  -L^z)  + 

must  vanish  at  upstream  and  downstream,  infinity’jj 

have 


so  that  we  must 


%  i/-  (■5/2)'^%  +X  =  9fC(l^'-Ui)(U~U^) 


(this  identity  can  of  course  be  derived  directly  from  the  shock 
conditions  IIIc9) 

From  this  identity  we  also  get 
{Z/Z)i'ilWi) 

and  the  equation  for  the  shock  structure  then  reduces  to 


^  n[it-  a/5)(Ui  hUz)]in-  o/^duitu^)] 


(III. .11) 


It  is  important  to  note  that  we  havej)  in  view  of  IIIplO  and 
assuming  to  represent  the  upstream  supersonic  conditions  ; 

OX-f 

and  therefore 

(2/5)CW,^U.2)>TA,  >U2  >  (l/o)(w,+u^)< 

so  that  the  denominator  of  III„11  is  always  positive  for 
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This  is  typical  for  3.  For  lesser  values  of  ^  ,  as 
always  encountered  in  physical  situations,  this  condition  does  not 
hold,  in  general,  and  one  can  expect  a  sufficiently  strong  wave  still 
to  exhibit  purely  kinematic  discontinuity  in  its  fine  structure  in 
the  absence  of  viscosity,  as  shown  in  ref.  13.  In  the  present  case, 
however,  the  velocity  profile  is  continuous  and  given,  implicitly 
by  simple  integration  of  III. 11,  i.e.  (using  III. 10) 

-Y+  I  In  (w,  -'i^)  —  Const.  (Ill .  12 ) 

The  first  member  being  a  monotonic  decreasing  function  of  TL  ,  the  plot 
of  U/  as  a  function  of  x  can  easily  be  derived  and  we  have  U 
for  X  =-fodand  U  =  U^f or  x  =  -oo  ,  which  indicates  the  irreversibility 
of  the  shock-process. 


APPENDIX  IV  -  Application  of  the  Fourier  Transformation  to  the  Solution 
of  the  Spatially  Homogeneous  Case  : 


We  use  the  transformed  equation  (7.15),  In  the  spatially 
homogeneous  case,  the  derivative  with  respect  to  %  disappears.  Further¬ 
more  we  can  always  take  0  and  hence,  =0,  because  %Lq  is  a 

constant  in  the  present  case  and  may  therefore  be  reduced  to  zero  by 
means  of  a  Galilean  transformation  which  does  not  change  the  spatially 
homogeneous  character  of  the  problem.  The  term  , which 

arises  because  of  the  linearization  does  not  occur  here  because  the 
spatially  homogeneous  case  is  essentially  linear. 


The  equation  to  be  treated  then  reduces  to  : 


(IV.  1) 


We  first  solve  the  problem  using  the  eigenvalue  method,  which 
reduces  here  to  the  computation  of  the  relaxation  times  and  the 

corresponding  relaxation  modes.  We  thus  have  to  look  for  solutions  of 
(IV. 1)  of  the  type 

(IV. 1)  then  reduces  to  an  ordinary  equation  for  > 


6J  4^  -f  =  o 


(IV. 2) 


This  represents  an  eigenvalue  problem  in  the  parameter  CT, 
the  corresponding  eigenfunction,  as  will  now  be  shown. 


representing 


The  general  solution  of  the  first  order  equation  IV. 2, 
member  is  immediately  found  to  be 


without  second 


(IV. 3) 


A  being  a  constant. 
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The  eigenvalue  nature  of  the  problem  results  from  the 
condition  that  will  be  the  Fourier  transform  of  a  function 
which  is  expected  to  behave  like  0  ^  for  large  1^  1  .  Therefore 

fiw)  exists  for  all  finite  complex  values  of  Cd  >  because 

{(Q  dc 

exists  even  if  cO  has  an  imaginary  part  .ytCj)is  thus  seen  to  be  an 
integral  analytic  function  and  this  is  a  strong  condition  on  the 
solution  IV. 3.  In  fact,  this  solution,  in  general  exhibits  a  singularity 
at  C0=  0  except  for  zero  or  negative  integral  values  of  the  parameter  Cr 
This  remark  solves  our  eigenvalue  problem,  for  we  must  have 


CT  =  for  0,  1,  2,  .  . 

The  corresponding  eigenfunctions  are 


\j/  =  const  Cu  S  (V  =<?'/, 2^...)  (IV. 4) 


Note  that  the  auxiliary  conditions 

/7x<iZ=o  / ~o 

translate,  after  Fourier  transformation,  into  the  simple  conditions 

/or  0J*o  J,  for  co=o 

These  conditions  are  obviously  satisfied  for  V  3,  so  that  the 
values  0,1  and  2  must  be  rejected  in  IV. 4, 

As  shown  in  ref  12,  the  inverse  Fourier  transforms  of  these 
eigenfunctions  are  the  Hermite  functions  indicated  in  paragraph  6. 

We  now  give  an  alternative  treatment,  using  Laplace 
transformations  with  respect  to  Z  and  based  on  a  mathematical  technique 
of  interest  for  more  complicated  cases-The  final  formula  so  obtained 
directly  shows  the  positive  definiteness  of  the  distribution  function. 
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If- we  apply  the  Laplace  transformation  with  respect  to 

f 

OJ)  dZ 

to  the  function  (p  , equation  LV.l  becomes 

jo  (p  ^  (pc O,0J)  -h  CO  ^ 


(IV. 5) 


Again,  the  general  solution  of  this  ordinary  differential 
equation  in  (0  is  immediately  obtained  as 

/  —CO /Z  (  1  — P  ”P  /  ^ b-l.r,) 

(p^e  j  Aoj  -hCJ  J  (p{o,CA^)e 


The  lower  integration  limit,  which  is  arbitrary,  has  been 
taken  as  zero  for  convenience  (as  will  appear)  and  this,  of  course, 
corresponds  to  a  particular  choice  for  the  constant  A. 

Just  as  in  the  eigenvalue  problem  treated  above, ^  must  be  an 
integral  analytic  function  of  the  complex  variable  6J,for  all  values 
of  p.  This  condition  will  define  the  constant  A, 

^£^60)  being  an  integral  analytic  function  of  6J,  we  may  write. 


cpCo,oS)  B 


^/z 


CO 


Z 


V-o. 


(IV, 7) 


The  expression  being  valid  for  all  Cd  , 


iV.7,  then  becomes 

=  e 


to  ^  ~h^ 


VrC* 


J 


(IV. 8) 


therefore 


For  arbitrary  p,  the  only  singularity  appears  in  the  A  term, and 

(1)  ^ 


A  =  0 


(1)  The  simplicity  of  this  condition  of  course  results  from  the  special 
choice  of  lower  limit  of  integration  in  formula  iV.6. 
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Going  back  to  the  original  function  (p  ,  we  get,  using  IV. 7 

co^{e'-i) 

Cp  ^  p Coe  ^  (IV. 9) 


The  inverse  Fourier  transformation  then  yields  the  following 
integral  representation  for  j-  T)  • 

icjz:-o-e-‘VO/z 


p(o^  cje'V 


■i-na  ^ 

..2n.  ^ 


^  r  Cojt -0-e''Vcoyz 

=  i/  ^ 

-Po 


rje" 


the  latter  expression  resulting  from  the  definition  of  (p  (0,00)  as 
Fourier  transform  of  ^(0,^). 


The  above  formula  may  be  written  as 


(IV. 10) 


~o0 


with  the  kernel 


Ccjfr-ze~^)—  — V^- 


-00 


which  is  the  Green  function  for  the  problem. 


Introducing  the  variable 


=  CO  aJ  i  —  Q 


-2C 


the  kernel  can  be  written 


\<(xtX= 


2(1-6-’-^' 


~t\2. 


e  —1 


i-e 


-zv 


2i\f 


00 


e'"^)  (IV.  11) 


is  clearly  a  positive  def initekernel  and  IV.  10  thus  shows 
that  ^  will  always  be  positive  if  it  was  initially  positive  for 
all  values  of  ^ 
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